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We have observed the cosmic microwave background temperature fluctuations 
+-> ■ in eight fields covering three separated areas of sky with the Very Small Array at 

34 GHz. A total area of 101 square degrees has been imaged, with sensitivity on angular 
scales 3?6-0?4 (equivalent to angular multipoles £=150-900). We describe the field 
selection and observing strategy for these observations. In the full- resolution images 
(with synthesised beam of FWHM ~ 17 arcmin) the thermal noise is typically 45 /iK 
and the CMB signal typically 55 ^k. The noise levels in each field agree well with the 
expected thermal noise level of the telescope, and there is no evidence of any residual 
systematic features. The same CMB features are detected in separate, overlapping 
observations. Discrete radio sources have been detected using a separate 15 GHz survey 
and their effects removed using pointed follow-up observations at 34 GHz. We estimate 
that the residual confusion noise due to unsubtracted radio sources is less than 14 
mJy beam -1 (15 fiK in the full-resolution images), which added in quadrature to the 
thermal noise increases the noise level by 6 %. We estimate that the rms contribution 
to the images from diffuse Galactic emission is less than 6 fj,K. We also present images 
which are convolved to maximise the signal-to-noise of the CMB features and are 
co-added in overlapping areas, in which the signal-to-noise of some individual CMB 
features exceeds 8. 
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ABSTRACT 



1 INTRODUCTION 

The measurement of primordial structure in the cosmic mi- 
crowave background (CMB) is of great importance in cos- 
mology. These features, imprinted at z ~ 1000, enable direct 
investi gation of the formation of structure in the universe as 



well as constraining the values of basic cosmological quanti- 
ties such as the amounts of various forms of matter and vac- 
uum energy in the universe. Primordial CMB fluctuations, 



however, are extremely faint (AT < 100/^K) and have to be 
seen through various foregrounds and in the presence of the 
systematic errors inevitably present in any experiment. 

These challenges account for the gap of nearly three 



decades b etween the discovery of the CMB ( Penzias fc Wil, 
sor [l965| ) and th e detection of prim ordial fluctuations by the 



COBE satellite flSmoot et al.||l992[ ). COBE detected statisti- 
cal anisotropy on angular scales of 90° to 10° (corresponding 
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Figure 1. Left: The array design used for the compact array. Right: The resulting Mu-coverage for a 6-hour obser- 
vation. 



to angular multipoles i = 2 to 20); this is greater than the 
horizon scale at the epoch of imprinting and provided very 
strong evidence for a period of inflation in the early uni- 
verse. After COBE, observations have focused increasingly 
on the angular scales at which features due to ac oustic os- 



cillat i ons in the primordial plasma were expected ( gakharov 
aunyaev 'doviclj "~ 



1965 



197C 



Peebles & Yu 



1970) 



Recently, results have been announced from a new 
generation of instruments capable of detecting a range of 
acoustic modes. As well as unequivocally detecting the firs t 
acoustic peak at £ 



220, MAXIMA (Lee et al 



20011 



BOOME RanG ( |Netterfield et alj [200 lj ) and DASI flHalver- 
son et |I] [2002j ), all detect power at smaller angular scales 
300 < I < 800, wit h strong evidence for a second peak at 
I ~ 550, while CBI flPadin et al||200l| ) (which is optimised 
for observation at high tj detects the sharp d ecline in p ower 
at I ~ 1000 expected from photon diffusion (Silk 1968). 



Here we report measurements of the CMB anisotropies 
on scales of 3? 6-0? 4 (angular multipoles ^=150-900) at 
34 GHz using the Very Small Array (VSA). This paper is 
the second in a series of four papers which report the re- 
sults of the first season of observations made using the VSA 
in its compact configuration. Here we describe the obser- 
vation al strategy, dat a redu ction and image-plane analysis, 
(2002) (hereafter Paper I) provide a 



Watson et al 



whilst 

detailed description of the experimental method and design. 
Extraction of the angular power spectrum is presented by 
Scott et al. (2002) (Paper III), and estimates of th e cosmo- 



logical param e ters u sing the VSA data are given by Etubino- 
Martin 



et al 



(|2002|) (Paper IV). A detailed description of 
the telescope will be given by Rusholme et al. (in prep.) 
(Paper 0). 



2 OBSERVATIONS 

The VSA is a development o f the Cosmic Anisotropy Tele- 
scope (CAT) (Robson et al 1993), a three-element inter- 
ferometer which operated at 15 GHz from a sea- level site 
in Cambridge. The VSA has 14 horn-reflector antennas 
mounted on a tilt-table and operates over the frequency 
range 26-36 GHz with an observing bandwidth of 1.5 GHz. 
The telescope is sited at the Teide Observatory, Tenerife at 
an altitude of 2400m. 

Prior to our first season of CMB observations, a series 
of commissioning observations was undertaken (see Paper 
I). These demonstrated that the telescope was working to 
specification and that systematic effects such as the effects 
of the Sun and the Moon can all be removed from the data 
to a very low level. Calibration of the data was shown to 
be limited only by the uncertainty in the absolute flux mea- 
surement of our primary calibrator, Jupiter, which is 3.5 %; 
the absolute flux c alibration of VSA o bservations is based 
on the flux scale of 



Mason et al. 



( 1999 ) and the brightness 



temperature of Jupiter at 34 GHz is taken to be 154.5 K. 

Detailed descriptions of the VSA, its data analysis 
procedures and commissioning observations are given in 
Rusholme et al. (in prep.) and Paper I; in this paper we 
describe the key features of the first season's observation of 
the CMB. 



2.1 Array configuration 

The 14 VSA antennas can be placed anywhere on the tilt- 
table, allowing freedom to design the array configuration 
for specific observational goals. The 'compact array' con- 
figuration used for the observations reported in this paper 
was optimised to provide almost uniform uu-coverage over 
the multipole range I ~ 150 - 900, whilst minimising the 
number of visibilities with low fringe rates. The sky signal 
in such visibilities cannot be separated from the spurious 
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Table 1. VSA field positions and total effective integration time remaining 
after flagging and filtering of the data. 





RA (J2000) 


DEC (J2000) 


Total Effective Integration Time (hrs) 


VSA1 


00 22 37.3 


30 16 38 


193 


VSA1A 


00 09 55.8 


30 33 12 


226 


VSA1B 


00 15 25.9 


28 04 40 


68 


VSA2 


09 37 57.0 


30 41 29 


271 


VSA2-OFF 


09 47 15.2 


30 41 05 


226 


VSA3 


15 39 39.3 


44 50 21 


187 


VSA3A 


15 30 35.1 


42 37 49 


206 


VSA3B 


15 46 20.8 


42 22 15 


119 



signal (see Paper I), and is predominantly contributed by 
the shortest north-south baselines. Figure 111 shows the VSA 
'compact array' configuration and uu-coverage. 



2.2 Observing frequency 

The VSA has the flexibility to observe anywhere in the 26- 
36 GHz range, with 1.5 GHz instantaneous bandwidth. This 
tunability allows the VSA to observe at more than one fre- 
quency. In principle, this allows one to fit for a component 
of a known spectral index, i.e. it would allow a separation 
of one or more Galactic components by virtue of their spec- 
tral indices. We have in fact chosen to observe at just a 
single frequency of 34 GHz. Our knowledge of the CMB 
foregrounds suggests that the Galactic foregrounds at fre- 
quencies ~ 30 GH z wi ll contribute a few fiK of signal in our 
fields (see section 3.2). Choosing a frequency at the higher 
end of the VSA observing band (34 GHz) reduces the free- 
free and synchrotron foregrounds by a factor of 1.8 — 2.2 
compared to the lower end (~ 26 GHz). Since the proposed 
s pinning dust component is expected at ~ 15 — 20 GHz 
(Draine & Lazarian 1998), its 34-GHz emission also will be 
considerably lower than its peak value. 



2.3 Field selection 

For practical CMB observation, it is important to choose 
fields which are relatively free from Galactic and extragalac- 
tic foregrounds. All the VSA fields are situated at Galactic 
latitudes greater than 20° and have low Galactic synchrotron 
and free-free emi ssion, as predic t ed by the 408-MHz all-sky 
radio survey o f 



Haslam et al 



(1982). The dust maps of 



Schlegel et al. (1998) were used to select fields with rela- 



tively low dust contamination. The actual level of Ga lactic 



contamination in our fields is discussed in Section 3.2 

To avoid bright clusters, we consulted the NO- 



195c 



RAS (IBohringer et alJ |200C|), BCS ( |Ebeling et al.|[1998| ) and 
Abell ( |ATelTn ~ 



catalogues. All fields were also chosen 
to be as free as possible of bright radio sources, since these 
are the major contaminant of CMB observations at 34 GHz. 



We used two low- frequency surveys , NVSS ( [Condo n et al 



199c) at 1.4 GHz and Green Bank (Gregory et al.'l996) at 



4.85 GHz to select CMB fields in which there were predicted 
to be no sources brighter than 500 mjy at 34 GHz within 
the VSA primary beam (FWHM= 4? 6). Predictions were 
made by extrapolating the flux density of every source in 
the 4.8 GHz catalogue to 34 GHz using its spectral index 



between 1.4 and 4.85 GHz. A further practical considera- 
tion which affected the choice of CMB fields was the need to 
observe all fields for a reasonable length of time (more than 
5 hours) from both Tenerife and Cambridge. This limits the 
declination range of our fields to between +26° and +54°. 
To maximize observing efficiency, we also selected fields that 
are evenly spaced around the sky. In order to increase the i- 
resolution of our measurement of the CMB power spectrum, 
we selected regions of sky where we can mosaic several VSA 
pointings. In each region of sky, mosaiced fields are sepa- 
rated by 2? 75. Our final choice of fields used during the first 
year of observation is shown in Figure ^ and in Table [l]. 



2.4 Observing strategy 

During the first year of VSA observations, we undertook two 
distinct observation programmes, each using the compact ar- 
ray. First, we have made deep mosaiced observations of eight 
fields in three evenly spaced regions of sky (hatched regions 
in Figure []). Each mosaiced field was observed for ~ 400 
hours, reaching a thermal noise of approximately 30 mjy. 
Mosaicing in this way enables us to increase the ^-resolution 
of our measurements whilst also reducing sample variance. 
However, all the fields must be observed with the Ryle Tele- 
scope in Cambridge for source subtraction (see section 3.5), 
and the time taken for us to survey all fields with the Ryle 
Telescope prior to any observation with the VSA has limited 
the area of sky that we can cover with deep mosaicing in this 
first year. Consequently, for I ^ 300, our measurements of 
the CMB power spectrum are limited by sample variance. 

In order to achieve a good estimate of the CMB power 
spectrum in the region i ^ 300, we undertook a second ob- 
servation programme; a shallow survey in one large area of 
sky. The shallow survey consists of 2 days observation on 
each of 30 mosaiced fields (Figure ^ hexagonal region) , and 
covers an area of approximately 300 square degrees. For this 
shallow survey, our source subtraction strategy no longer 
limits the area of sky we can observe with the VSA in a 
given time. Since for the shallow survey we are only con- 
cerned with low-£ observations, where the contribution of 
point sources to the CMB power spectrum will be negligible, 
prior surveying with the Ryle Telescope is not necessary. In- 
stead, we monitor only sources predicted to be greater than 
100 mjy at 34 GHz on the basis of low-frequency survey 
information. Only the results from the deep mosaiced fields 
are reported in this series of VSA results papers; the results 
from the shallow survey will be presented in a later paper. 
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Figure 2. Plot of the VSA fields observed during our first year of 
observations, projected on to the equatorial plane. The hatched 
regions are fields which were observed using deep mosaicing, the 
results of which are the subject of this series of papers. The re- 
maining fields were observed as part of a shallow survey and are 
not discussed further here. Radio sources predicted to be brighter 
than 500 mjy at 34 GHz are displayed as star symbols. A cut-off 
at Galactic latitude |b| > 20° has been applied. 



which account for abnormal errors such as pointing, hard- 
ware failure, bad weather and filtering of the spurious signal. 

The fully-flagged and stacked data are held as visibility 
files, containing the real and imaginary part for each ob- 
served uu-position along with an associated rms noise level. 
It is critical for the power spectrum analysis that the es- 
timate of the noise on each visibility is accurate. We cal- 
culate the noise on each visibility point from the scatter 
in the visibility data on each baseline each day. This takes 
into account the number of flagged points in each smoothed 
64-second sample, in addition to other flagging operations 
such as gain corrections implied by the system temperature 
monitoring. As a further check, we recalculate our noise es- 
timates when the data are gridded in the w-plane, prior to 
producing reduced datasets for the power spectrum estima- 
tion. Here we estimate the noise for each gridded visibility 
using the spread of values within the stacked files (i.e. using 
the complete 70-80 day data sets). Both methods yield the 
same noise estimates. 

The final data stacks for each of the VSA fields typically 
contain ~ 800,000 visibility points, each one averaged over 
64 seconds. The combined edits from all categories result in 
a rejection of ~ 50% of the data in each field. 



2.5 Observing schedule 

Quasi-independent tracking of the VSA antennas enables 
us to effectively filter out the effects of bright contaminat- 
ing sources such as the Sun and the Moon (see Paper I). 
We are able therefore to observe during both the day and 
night, with useful CMB observations possible at separations 
as close as 40 degrees from the Sun and 30 degrees from the 
Moon. Consequently, we implemented an observing strategy 
that made use of all available time during the first year of 
observation. 

In each 24-hour period, we obtained 5-6 hours of in- 
tegration on each of 3 CMB fields; the total effective in- 
tegration time (after flagging and filtering of the data) for 
each of the VSA fields is given in Table bL The remaining 
time in each 24-hour period was divided between calibration 
sources. We observed at least one flux calibrator (Jupiter, 
Cas A or Tau A), with observations of fainter sources (e.g. 
NGC7027, 3C273 and Cyg A) interspersed. Observation of 
fainter sources allowed valuable checks on the diurnal sta- 
bility of the telescope. 



2.6 Data reduction 

The data from the first year of VSA observations were re- 
duced as described in detail in Paper I. Each day's CMB 
observation is reduced and analysed independently. The re- 
sulting calibrated visibilities from each day are stacked to- 
gether, and the full dataset for each field is analysed. We 
check the stacks by eye for signs of residual contaminating 
signal, which is easily identified as regions of coherent phase. 
The data are stacked both as a function of hour angle and 
per-antenna, and are nagged interactively. This stage of flag- 
ging typically removes 20% of the data; this is in addition 
to the preliminary flagging operations described in Paper I, 



3 FOREGROUNDS 

3.1 Diffuse Galactic foregrounds 

The diffuse Galactic foregrounds consist of three known 
components: synchrotron, free-free and vibrational dust, to- 
gether with a further postulated component of spinning 
dust. 

Synchrotron emission ha s a varying spectral index (T oc 
j/ -/3 ) of (3 aynch 2.7 — 3.2 (Davies et al. 1996] ) depending 
on frequency and position on the sky. The steep spectral 
index allows low-frequen cy radio surveys su ch as the 408- 
MHz all-sky rad io map ([Haslam et al.| |l982j ) and the 1420- 
MHz radio map (Reich & Reich 1988) to be used to estimate 
the synchrotron component at higher frequencies. 

Free-free emission is thermal bremsstrahlung radiation 
emitted by ionized gas in the Galaxy. The spectral index 
for free-free is more well-known with /3ff ~ 2.15; it shows 
very little variation with electron temperature and density, 
although it does steepen slightly at higher frequencie s. The 
656.3 nm Ha line is a good tracer of free-free emission ( Valls- 



Gabaud ( 199S ) and Dickinson et al. in prep). Dust absorp- 
tion, however, can cause these estimates to be systematically 
lower than the true values. 

Vibrational dust is only a concern for C MB experiments 
obser ving at frequencies above ~ 90 GHz (ie Zotti et al 



199S), and is expected to contribute << 1 fiK at VSA ob- 



serving frequencies. A further component has recently been 



postulated by Draine & Lazarian ( 199S ) , based on observa- 


tional evidence (Kogut et al. 


(1996); Leitch et al. (1997); 


de Oliveira-Cost 


a et al. (1997, 1999, 2005), but see also 


Mukherjee et al 


(2001)). They postulate a population of 



rotational emission. The spectrum of spinning dust is likely 
to be highly peaked at about 15 — 20 GHz. 
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Table 2. Galactic foreground contamination for the 3 VSA mosaiced fields 
at 1° resolution. The rms estimates are in fiK, and are calculated for 34 GHz. 



Component 


Template / assumptions 


VSA1 


VSA2 


VSA3 


Synchrotron 


Haslam 408 MHz 


2.2 


1.2 


1.0 




P = 3.0 








Free-free 


WHAM Ha 


0.8 


0.8 


0.6 




1R=4.5 /jK 








Spinning Dust 


Schlegel et al. 100 lira 


5.0 


1.4 


2.4 




1 MJy sr" 1 = 10 








Total (estimated) 


uncorrelated 


5.5 


2.0 


2.7 



3.2 Galactic foregrounds in the VSA fields 

A simple method of estimating the foreground contamina- 
tion is to calculate the rms variation in each of the fore- 
ground template maps for each VSA field. This is then con- 
verted to temperature units at 34 GHz. This is simplistic 
and susceptible to uncertainty due to the extrapolation to 
34 GHz. This method, however, does give us an estimate 
of the amplitude of the foreground signals and indicates to 
what level the VSA observations are contaminated by diffuse 
foregrounds. A more rigorous cross-correlation analysis will 
be presented in a future paper (Dickinson et al., in prep.). 

The rms method has been applied to 8° x 8° patches 
of sky centred on each of the three VSA deep fields (VSA1, 
VSA2 and VSA3). These are large enough to cover the mo- 
saiced fields within the 4? 6 FWHM of the primary beam. 
We use the de-striped and source-subtracted version of the 
Hasla m et al. (1982) 408 MHz map as used by 



Davies et al 



(1996) to trace the synchrotron com ponent, whilst the re - 



cently published WHAM Ha data (Reynolds et al 



199£ 



are used to trace the free- free component. For the dust- 
correlate d emission, we use th e 100 fim IRAS/DIRBE map 
given by schlegel et al. (1998). The maps were smoothed to 
1° resolution since the WHAM data are at ~ 1° resolution. 

The synchrotron emission is assumed to have a spec- 
tral index (3 = 3.0 between 408 MHz and 34 GHz. The 
conversion from Hcv units (1 Rayleigh (R) = 10 6 /4-7r pho- 
tons s" 1 cm" 2 sr" 1 ) to free-free emission is calculated to 
be 4.5 fiK/R (Dickinson et al. in prep.). To estimate the 
total dust-correlated component (which may include emis- 
sion from spinning dust and vibrational dust) we have as- 
sumed a coupling coefficient of 10 /^K/(MJy sr -1 ). This 
number is rather arbitrary since previous estimates which 
range from — 36/iK/(MJy sr" 1 ) at ps 30 GHz, have been 
calculated for different angular scales and for different re- 
gions of sky. The total value is the sum of the individual 
components added in quadrature. The estimates of the com- 
ponents and total contamination from Galactic foregrounds 
for the three fields are given in Table ^. 

The value of the synchrotron component is probably an 
over-estimate since there are some residuals from discrete 
radio sources left in the map; above ~ 10 GHz, the free-free 
component is expected to dominate the synchrotron compo- 
nent. The free-free estimate from Ha is probably the most 
reliable although there may be some absorption of the Ha 
signal due to dust. Following Dickinson et al. (in prep.), we 
can use the 100 /im emissivity to calculate what this correc- 
tion might be. For all the VSA fields, the dust emissivity is 
< 5 MJy sr" 1 . If all the dust were lying in front of the Ha 



gas, then the correction is 23 percent. A more likely situa- 
tion is a mixing of the gas and dust, thus the dust correction 
is ~ 10 percent. The table indicates the brightest Galactic 
foreground may be the spinning dust, although very little is 
known about this component. The estimates are likely to be 
higher than the real value since we chose a rather high cou- 
pling coefficient based primari l y on the COBE-DMR dust- 
correlated value (Kogut et al, 1996) which is only strictly 



true at the 7° scale. This is the total dust correlated compo- 
nent which will also include any other correlated foreground, 
most notably, the free-free component. 

In summary, using rms estimates for these VSA fields, 
it seems that the diffuse Galactic foregrounds are very weak 
at 34 GHz compared to the primordial CMB signal. For 
example, if Tomb — 60 fiK near the first peak I ~ 200 and 
Tgal = 5 fiK, then the contamination is at the 0.35 per cent 
level. For t ~ 800, if Tomb = 30 /iK, then the contamination 
is 1.4 per cent. However, the contamination is likely to be 
considerably smaller at higher ^-values, since all the Galactic 
foregrounds have power spectra which decrease at smaller 
angular scales. The full cross-correlation analysis (Dickinson 
et al., in prep.) will provide a more accurate estimate of the 
contamination from Galactic foregrounds. 



3.3 Point sources 

Radio sources (mostly radio galaxies and quasars, but also 
some stars) are expected to be a major contamina nt of 



CMB obs e rvatio ns at centimetre wavelengths (see e.g., Tay 
(2001)). Their population is not well known at 



lor et al 



frequencies higher than 10 GHz, and many are expected to 
be variable and/or have spectra that rise with frequency, 
i.e. a < (where flux density S oc v~ a ). On the angular 
scales of interest to primordial CMB work, such sources are 
also generally unresolved. To eliminate point source con- 
tamination from VSA observations, all the radio sources in 
each CMB field with flux density greater than our source- 
confusion limit, are observed and subtracted directly from 
our visibility data. We estimate the contribution of unsub- 
tracted sources below our confusion limit using a prelimi- 
nary estimate of the 34 GHz source count; this source count 
is estimated from the sample of sources we detect at 34 GHz. 



3.4 Source confusion and the VSA 

The level to which sources must be subtracted from our 
CMB fields is determined by comparing the flux sensitivity 
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Table 3. Point sources subtracted from VSA fields. The apparent 
flux is the flux of each source at 34 GHz as observed by the VSA 
(i.e. including the effect of the VSA primary beam). This flux 
represents the actual flux subtracted from each map. 

RA (J2000) DEC (J2000) Apparent Flux (Jy) 



Table 3 - continued 



RA (J2000) DEC (J2000) Apparent Flux (Jy) 



VSA1 


00 


10 


27.7 


28 


54 


58 


0.028 




00 


15 


06.1 


32 


16 


13 


0.155 




00 


19 


39.7 


26 


02 


53 


0.043 




00 


23 


39.6 


29 


28 


29 


0.091 




00 


24 


34.3 


29 


11 


31 


0.088 




00 


28 


17.0 


29 


14 


29 


0.054 




00 


34 


43.4 


27 


54 


26 


0.037 




00 


36 


59.7 


29 


58 


59 


0.029 




00 


41 


15.9 


32 


11 


08 


0.011 


VSA1A 


23 


55 


44.0 


30 


12 


42 


0.018 




00 


10 


27.7 


28 


54 


58 


0.058 




00 


12 


47.3 


3.3 


53 


36 


0.029 




00 


13 


44.0 


34 


41 


42 


0.008 




00 


15 


06.1 


32 


16 


13 


0.167 




00 


15 


37.9 


33 


32 


11 


0.017 




00 


19 


37.7 


29 


56 


02 


0.033 




00 


23 


39.6 


29 


28 


29 


0.019 




00 


24 


34.3 


29 


11 


31 


0.020 


VSA1B 


00 


00 


35.1 


29 


14 


29 


0.019 




00 


00 


55.9 


25 


16 


24 


0.014 




00 


01 


21.5 


25 


26 


51 


0.015 




00 


0.3 


16.3 


24 


59 


42 


0.008 




00 


04 


34.5 


29 


46 


18 


0.021 




00 


06 


48.8 


24 


22 


38 


0.027 




00 


12 


38.1 


27 


02 


10 


0.057 




00 


15 


06.1 


32 


16 


13 


0.033 




00 


16 


03.6 


24 


40 


18 


0.025 




00 


18 


12.4 


29 


21 


25 


0.101 




00 


19 


39.7 


26 


02 


5.3 


0.298 




00 


19 


52.4 


26 


47 


31 


0.052 




00 


2.3 


23.0 


25 


39 


19 


0.044 




00 


34 


43.4 


27 


54 


26 


0.015 


VSA2 


09 


2.3 


48.0 


31 


07 


56 


0.047 




09 


23 


51.5 


28 


15 


26 


0.047 




09 


25 


43.6 


31 


27 


11 


0.062 




09 


32 


55.0 


3.3 


39 


29 


0.023 




09 


39 


01.6 


29 


08 


29 


0.055 




09 


11 


48.1 


27 


28 


38 


0.051 




09 


52 


06.1 


28 


28 


31 


0.011 




09 


53 


27.9 


32 


25 


52 


0.012 




09 


58 


20.9 


32 


24 


02 


0.014 



of the VSA to the rms confusion noise remaining from un- 
subtracted sources. The VSA, operating in its compact con- 
figuration, has a flux sensitivity of about 30 mjy in each deep 
mosaic pointing. In order that our observations are domi- 
nated by receiver noise we must therefore ensure that the 
rms flux density on a VSA map from unsubtracted sources, 
the confusion noise, is less than the flux sensitivi t y of 3 mjy. 
Following the standard approach (e.g. Scheuer (1957)), 



dconf = fi / n(S)S 2 dS, 
Jo 



where S is the flux density, n(S) is the differential source 
count per steradian, f2 is the VSA synthesised beam solid 
angle and S S ub is the flux density limit down to which sources 
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are subtracted. Normalising this in terms of the number of 
sources, N, per VSA synthesised beam that are subtracted, 

N= D.n{S)dS, 

and parameterising n(S) as n(S)dS = KS -13 , we obtain 

<W = [7V(/3-l)/(3-/3)] 1/2 S sub . 

For a given level of source subtraction, Ssub, one can 
thus estimate the residual confusion noise. The level of 
source subtraction can then be optimised to ensure that 
o"conf is less than the flux sensitivity of the VSA. To calculate 
this confusion noise, however, the differential source count, 
n(S), at our observing frequency of 34 GHz is required. At 
the outset of observations with the VSA, the source count 
at 34 GHz was not known. Instead we used an estimate 
extrapolated from the 15-GHz source count obtained using 
source subtraction observations for the Cosmic Anisotropy 
Telescope(CAT). 
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Table 4. This table gives the thermal noise level, in mjy beam -1 , for each VSA field, measured both from standard 
and autosubtracted maps (see text for explanation) . The rms noise level in the centre of each map and the corre- 
sponding residual CMB signal observed in each field, in mjy, is also given. These values are taken from CMB maps 
made at full resolution. The approximate flux-to-temperature conversion is 1 mjy beam -1 = 1.1 fiK, giving mean 
rms CMB fluctuations in the maps of 40 - 65 fiK. 
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Based on the small sample of sources observed during 
the CAT source-subtraction progra m, the 15 GHz s ource 
counts were found to have j3 = —1.8 ( |0' Sullivan 1995). Us- 
ing this value for /3 at 15 GHz and assuming a conservative 
average spectral index of a~0 between 15 and 34 GHz, we 
find that to ensure that the confusion noise is significantly 
lower than the flux sensitivity of the VSA (~ 30 mjy) we 
must subtract all sources brighter than 80 mjy at 34 GHz. 



3.5 Source-subtraction strategy 

To ensure that all potentially contaminating sources are de- 
tected, one would ideally survey each CMB field at the same 
frequency as the VSA, but with increased resolution and sen- 
sitivity. This, however, would require a telescope more ex- 
pensive than the VSA. Instead, we adopt a two-stage strat- 
egy- 
First, prior to observation with the VSA, we survey all 
the VSA fields at 15 GHz using the Ryle Telescope (RT) 
in Cambridge. The RT, which uses five 13-m diameter an- 
tennas and gives a resolution of ~ 30 arcsecs, is used in 
a raster scanning -mode and reach e s an rms noise level of 
a ~ 4 mjy (see Waldram et al. (2002)). This allows us 
to identify sources above 20 mjy at 15 GHz, and ensures 
that we find all sources above our source-confusion limit of 
80 mjy at 34 GHz, even allowing for a spectral index as 
steep as a ~ —2 between 15 and 34 GHz. 

Having identified the contaminating sources in each 
field, we monitor each source at 34 GHz using a separate 
single-baseline interferometer working simultaneously with 
the VSA (see Paper I). This single-baseline interferometer 
consists of two 3.7-m dishes separated by 9 metres on a 
north-south baseline. Each dish is situated in an enclosure 
similar to that of the VSA, and is fitted with identical horn- 
reflector feeds and receivers as used on the main array. All 
sources identified in the 15 GHz survey are monitored. The 
monitoring is done simultaneously with and at the same fre- 
quency as the VSA observation, ensuring that sources which 
are variable on time-scales as short as a few days can be sub- 
tracted accurately. 

On average, each source identified in the 15-GHz survey 
is observed once every three days. The flux of each monitored 
source is then averaged over the duration of each correspond- 



ing field observation (typically 80 days) and its contribution 
is then subtracted directly from the stacked visibility data. 



3.6 Point source observations 

In total, 346 sources with flux densities greater than 20 mjy 
were identified in the 15 GHz survey of the VSA fields. All 
of these sources were observed with the source subtractor, 
and the list of monitored point sources with flux densities 
greater than 60 mjy at 34 GHz is given in Table |S| Al- 
though our source-subtraction strategy requires that only 
sources brighter than 80 mjy must be subtracted, and our 
source list is only complete to this level, since the data were 
available we subtracted point sources down to a flux limit of 
60 mjy. Typically, between 10 and 15 sources with flux den- 
sity greater than 60 mjy were subtracted from each CMB 
field. 

To calculate the residual point-source contamination 
remaining after source-subtraction, we used a preliminary 
source count at 34 GHz. The source count is based on the 
78 sources with S34 > 60 mjy detected by the source sub- 
tractor. The differential count is given approximately by 



dN 
dS 



54 



T "I -1 

Jy sr . 



Note that this count is not entirely accurate due to the 
uncertain completeness at lower flux levels, but is sufficient 
to estimate our residual confusion noise. A more complete 
analysis of the 34-GHz count will be given in a subsequent 
paper. When the count is integrated up to our nominal sub- 
traction limit of 80 mjy it gives a residual confusion noise 
of (Tconf = 14 mjy, which when added in quadrature to our 
typical thermal noise of 40 mjy gives only a 6 % increase in 
noise level. The effect of the residual sources on the CMB 
power spectrum is discussed in Paper III. 



4 RESULTS 

Although the CMB power spectrum is calculated directly 
from the complex visibility data obtained from our observa- 
tions (see Paper III) , image-plane analysis provides valuable 
consistency checks on the data. Here we present images of 
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Figure 3. Maps of a sample VSA field (VSA3A) at full resolution. Top left: the map before source subtraction. 
Positions of sources to be subtracted are marked with crosses, and the FWHM of the primary beam is indicated by a 
circle. CMB fluctuations (and some sources) are clearly visible within the primary beam. Top right: the synthesised 
beam of all the images. Bottom left: the map after source subtraction. Bottom right: the autosubtracted map (see 
text), showing the uniform noise level across the map. In all maps the greyscale runs from -150 to 150 mjy beam -1 
and the contour interval is 50 mjy beam -1 ; in the beam the greyscale is -0.3 to 0.3 and the contour interval is 5% 
of the peak. 



the VSA fields and discuss the checks that we have per- 
formed on the data. 

Figure ^ shows images of a typical VSA individual field 
(VSA3A) at full resolution. The top two panels show the 
map before source-subtraction and the synthesised beam, 
while the lower left panel shows the map after source sub- 
traction. Inside the primary beam envelope (FWHM = 4? 6, 
shown by a circle), the maps are dominated by astronomi- 
cal signal. Outside this region, the maps are dominated by 
instrumental (thermal) noise. 

Comparison of the rms power measured in the outer 
region of each map with the calculated instrumental noise 



provides information about the level of residual contaminat- 
ing signals. Unwanted signals (such as crosstalk or distant 
bright sources) are not constrained to lie within the pri- 
mary beam envelope in the image plane. Instead their effect 
is to increase the total rms power level across the whole of 
a map. A simple way to measure the instrumental noise di- 
rectly from a map is to measure the rms level far from the 
primary beam. The measured noise levels, given in Table 
m all agree extremely well with the calculated values. The 
slightly higher noise levels obtained in VSA1B and VSA3B 
reflect the smaller amount of data collected on these fields. 

A more rigorous approach to estimating the thermal 
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noise on each map, and also the amount of any residual con- 
taminating signal, is to 'autosubtract' our visibility data. 
Here we take the time-ordered visibility data for each field 
and reverse the sign of alternate visibilities. On the time 
scale between adjacent visibility points (64 seconds), both 
astronomical and spurious signals are effectively coherent, 
and are thus cancelled out. In contrast, the noise on adja- 
cent visibilities is completely uncorrelated, and the rms noise 
on the map is unaffected. This technique therefore provides 
a robust estimate of our thermal noise. In addition, we can 
compare the rms power in the autosubtracted maps with 
that far outside the primary beam envelope of our standard 
field maps. Any discrepancy between the two would be in- 
dicative of residual spurious signals. 

Autosubtraction was applied to all eight VSA fields, and 
maps of each field were made. The bottom right panel in Fig- 
ure^ shows the autosubtracted map for VSA3A. Clearly the 
astronomical signal has been subtracted, leaving a constant 
noise level. The rms noise on each autosubtracted map was 
measured and the values are given in Table ^| There is excel- 
lent agreement between the noise levels measured far from 
the map centres and those measured by autosubtraction. 

The approximate CMB signal level in each map can be 
estimated by subtracting the map rms in the central area 
(4x4 degrees) in quadrature from the rms noise level esti- 
mated from the autosubtracted maps (or from the far-out 
portions of the sky maps). These figures are also given in 
Table ^, and represent the average temperature fluctuations 
in the CMB averaged over the Grange of the observations. 
For the synthesised beam size of the full resolution maps 
(~ 17 arcmin), the approximate flux-to-temperature con- 
version is lmjybeam -1 = 1.1 fiK, giving mean rms CMB 
fluctuations in the maps of 40 - 65 [iK. 

All the maps are robust to minor variations in flagging 
and filtering. The data have also been split by observing 
epoch, by day/night and by alternate visibilities; the result- 
ing maps are fully consistent with each other. Analogues 
of theses consistency checks are also applied to the power 
spectra (see Paper III). 

In addition to analysing each field independently, it 
is instructive to compare images of fields within each mo- 
saiced region. As an example, we present an image of the 
VSA2/VSA2-OFF mosaic (Figure Q). Here a contoured map 
of VSA2-OFF is overlaid on a grey-scale image of the VSA2 
field. Common structures are clearly seen, agreeing in both 
position and in flux density. 

Finally, in Figure ^| we present images of the fully 
source-subtracted VSA fields. To enhance signal-to-noise, 
a itu-taper function has been applied to the data. The ta- 
per has the form of the CMB power spectrum estimated 
from the complete VSA data set (see Paper III). The indi- 
vidual overlapping maps are combined linearly using stan- 
dard AlPS tasks. The individual maps are transformed to the 
same grid, then corrected for the individual primary beam 
responses (which multiplies up the outer parts of the map); 
they are then combined pixel by pixel, weighted by their pri- 
mary beams, and the resulting image multiplied by the total 
weight map, resulting in an image in units of signal-to-noise. 
This is then re-scaled back to the map units of Jy beam -1 . 
(This final re-scaling is only approximate as the sensitivi- 
ties vary slightly between the component maps; however, no 
quantitative analysis is done using these images). The peak 
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Figure 4. Comparison of the overlapping region in the VSA2 
mosaic. A contour map of VSA2-OFF is overlayed on a greyscale 
image of VSA2. Structures common in both position and flux level 
are visible. Both maps were made with a taper corresponding to 
the CMB power spectrum to maximise the signal-to-noise, have 
been corrected for the primary beam response of the VSA, and are 
truncated at the 25% point of the primary beam. The grey scale 
runs from —300 to 300 mjy beam -1 and the contour interval is 
150 mjy beam -1 . 



value of the signal-to-noise varies between 6 and 8.5 in the 
three combined images. The maps have not been CLEANed, 
so the CMB fluctuations are convolved with the synthesised 
beam, which is somewhat different for each field due to dif- 
ferences in flagging. 



5 CONCLUSIONS 

We have observed eight overlapping fields in three separate 
regions of the sky with the VSA compact array, a total of 101 
square degrees. The CMB anisotropies are clearly detected 
in all fields, and individual features in different overlapped 
pointings agree well with each other. We have assessed the 
Galactic foreground contamination based on low-frequency 
radio, dust, and H a data, and found it to be negligible. To 
eliminate confusion by foreground radiosources, we have sur- 
veyed all the fields at 15 GHz with the Ryle Telescope (RT). 
The flux of all sources found in the RT survey were moni- 
tored simultaneously with and at the same frequency as the 
VSA observations using a separate single-baseline interfer- 
ometer. These sources were subtracted from the data. The 
count of sources in the VSA fields is effectively complete 
to 80 mjy at 34 GHz, although we have detected sources as 
faint as 60 mjy. After removal of these sources from the VSA 
main array data, the residual source confusion is negligible. 
We have checked the data for evidence of components other 
than the CMB and thermal noise, and have found none. 
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Figure 5. Images of the three VSA regions, with radio sources subtracted. The individual maps were made with a 
radial taper applied to the visibilities in order to enhance the signal-to-noise of the CMB features, then combined 
using the AlPS tasks LTESS and STESS. These form a linear combination of the component maps, weighted by the 
individual primary beam responses. The point-spread function (or synthesised beam) have not been deconvolved 
from these images, which thus have correlations across them. Differing aperture plane coverages due to different 
flagging of the data have resulted in different synthesised beams in each image. Top left: VSA1,1A and IB. Top 
right: VSA2 and VSA2-OFF. Bottom left: VSA3.3A and 3B. In each case the greyscale runs from -300 to 300 mjy 

beam -1 . Bottom right: The weighting function used to taper the maps in order to improve the signal-to-noisc ratio, 

1/2 

based on the power spectrum of the data; the weighting function is proportional to C, . 



